## Loading required package: foreach
## 
## Attaching package: 'foreach'
## The following objects are masked from 'package:purrr':
## 
##     accumulate, when
## Loading required package: iterators
## 
## Attaching package: 'data.table'
## The following objects are masked from 'package:lubridate':
## 
##     hour, isoweek, mday, minute, month, quarter, second, wday, week,
##     yday, year
## The following objects are masked from 'package:dplyr':
## 
##     between, first, last
## The following object is masked from 'package:purrr':
## 
##     transpose
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.040006 -0.003320  0.000519  0.003040  0.120229 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0022478  0.0001405  -15.99   <2e-16 ***
## test_set$FCH4  0.1414673  0.0071422   19.81   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006864 on 2441 degrees of freedom
## Multiple R-squared:  0.1385, Adjusted R-squared:  0.1381 
## F-statistic: 392.3 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.155372 -0.002183  0.000236  0.001852  0.097756 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.299e-03  6.117e-05  -21.23   <2e-16 ***
## result$FCH4  4.859e-01  2.776e-03  175.03   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006003 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.7582, Adjusted R-squared:  0.7582 
## F-statistic: 3.063e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.046687 -0.003331  0.000832  0.003463  0.069299 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0029299  0.0001406  -20.84   <2e-16 ***
## test_set$FCH4  0.1077385  0.0069874   15.42   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006898 on 2441 degrees of freedom
## Multiple R-squared:  0.08875,    Adjusted R-squared:  0.08838 
## F-statistic: 237.7 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.170860 -0.002059  0.000387  0.002018  0.094157 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.614e-03  6.428e-05  -25.11   <2e-16 ***
## result$FCH4  4.712e-01  2.917e-03  161.52   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006308 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.7276, Adjusted R-squared:  0.7275 
## F-statistic: 2.609e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.049684 -0.003209  0.000549  0.003025  0.053550 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0023895  0.0001374  -17.39   <2e-16 ***
## test_set$FCH4  0.1334830  0.0067873   19.67   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006744 on 2441 degrees of freedom
## Multiple R-squared:  0.1368, Adjusted R-squared:  0.1364 
## F-statistic: 386.8 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.111600 -0.002135  0.000346  0.001931  0.091834 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.492e-03  6.094e-05  -24.48   <2e-16 ***
## result$FCH4  4.741e-01  2.765e-03  171.45   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00598 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.7506, Adjusted R-squared:  0.7505 
## F-statistic: 2.94e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.057383 -0.003116  0.000761  0.003353  0.091518 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0028336  0.0001358  -20.87   <2e-16 ***
## test_set$FCH4  0.0910863  0.0052743   17.27   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006696 on 2441 degrees of freedom
## Multiple R-squared:  0.1089, Adjusted R-squared:  0.1085 
## F-statistic: 298.3 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.166895 -0.002164  0.000562  0.002347  0.079910 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.938e-03  6.801e-05   -28.5   <2e-16 ***
## result$FCH4  4.025e-01  3.086e-03   130.4   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006674 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.6352, Adjusted R-squared:  0.6352 
## F-statistic: 1.701e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.050944 -0.003289  0.000543  0.003246  0.086515 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0023503  0.0001381  -17.02   <2e-16 ***
## test_set$FCH4  0.1255967  0.0072266   17.38   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006756 on 2441 degrees of freedom
## Multiple R-squared:  0.1101, Adjusted R-squared:  0.1098 
## F-statistic: 302.1 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.121750 -0.002120  0.000285  0.001945  0.103527 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.408e-03  6.142e-05  -22.92   <2e-16 ***
## result$FCH4  4.882e-01  2.787e-03  175.15   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006027 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.7585, Adjusted R-squared:  0.7584 
## F-statistic: 3.068e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.067230 -0.003142  0.000826  0.003442  0.056345 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0027802  0.0001475  -18.85   <2e-16 ***
## test_set$FCH4  0.0932542  0.0066440   14.04   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007261 on 2441 degrees of freedom
## Multiple R-squared:  0.07468,    Adjusted R-squared:  0.0743 
## F-statistic:   197 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.161268 -0.002152  0.000460  0.002190  0.085538 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.757e-03  6.767e-05  -25.97   <2e-16 ***
## result$FCH4  4.483e-01  3.071e-03  145.97   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006641 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.6857, Adjusted R-squared:  0.6856 
## F-statistic: 2.131e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.070367 -0.003234  0.000691  0.003299  0.040974 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0024532  0.0001389  -17.66   <2e-16 ***
## test_set$FCH4  0.1285094  0.0061965   20.74   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006815 on 2441 degrees of freedom
## Multiple R-squared:  0.1498, Adjusted R-squared:  0.1495 
## F-statistic: 430.1 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.169359 -0.002225  0.000399  0.002048  0.088535 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.487e-03  6.353e-05  -23.41   <2e-16 ***
## result$FCH4  4.518e-01  2.883e-03  156.69   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006235 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.7154, Adjusted R-squared:  0.7153 
## F-statistic: 2.455e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.041657 -0.003334  0.000564  0.003209  0.067046 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0025516  0.0001363  -18.73   <2e-16 ***
## test_set$FCH4  0.1411725  0.0058544   24.11   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006699 on 2441 degrees of freedom
## Multiple R-squared:  0.1924, Adjusted R-squared:  0.1921 
## F-statistic: 581.5 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.163063 -0.002234  0.000370  0.002103  0.060342 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.615e-03  6.224e-05  -25.95   <2e-16 ***
## result$FCH4  4.433e-01  2.824e-03  156.95   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006108 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.716,  Adjusted R-squared:  0.716 
## F-statistic: 2.463e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.058322 -0.003148  0.000593  0.003126  0.045491 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0025190  0.0001281  -19.67   <2e-16 ***
## test_set$FCH4  0.1084774  0.0061025   17.78   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006298 on 2441 degrees of freedom
## Multiple R-squared:  0.1146, Adjusted R-squared:  0.1142 
## F-statistic:   316 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.148409 -0.002143  0.000404  0.002087  0.078711 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.641e-03  6.185e-05  -26.53   <2e-16 ***
## result$FCH4  4.584e-01  2.807e-03  163.32   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006069 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.7319, Adjusted R-squared:  0.7319 
## F-statistic: 2.667e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.06683 -0.00382  0.00015  0.00262  0.38233 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0017515  0.0002697  -6.494 1.01e-10 ***
## test_set$FCH4  0.1675527  0.0128210  13.069  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01319 on 2441 degrees of freedom
## Multiple R-squared:  0.06539,    Adjusted R-squared:  0.06501 
## F-statistic: 170.8 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.19496 -0.00221  0.00009  0.00168  0.37161 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.092e-03  8.816e-05  -12.39   <2e-16 ***
## result$FCH4  5.268e-01  4.001e-03  131.67   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008652 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.6396, Adjusted R-squared:  0.6396 
## F-statistic: 1.734e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.052202 -0.003838  0.000688  0.003471  0.054760 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0028151  0.0001477  -19.07   <2e-16 ***
## test_set$FCH4  0.0859087  0.0071489   12.02   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007217 on 2441 degrees of freedom
## Multiple R-squared:  0.05586,    Adjusted R-squared:  0.05547 
## F-statistic: 144.4 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.166321 -0.002195  0.000274  0.001925  0.186303 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.362e-03  7.202e-05  -18.92   <2e-16 ***
## result$FCH4  4.980e-01  3.268e-03  152.39   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007067 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.7039, Adjusted R-squared:  0.7039 
## F-statistic: 2.322e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.037570 -0.003370  0.000515  0.003148  0.058573 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0025736  0.0001416  -18.17   <2e-16 ***
## test_set$FCH4  0.0835142  0.0066356   12.59   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006946 on 2441 degrees of freedom
## Multiple R-squared:  0.06094,    Adjusted R-squared:  0.06055 
## F-statistic: 158.4 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.191201 -0.002205  0.000319  0.001998  0.099815 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.519e-03  6.668e-05  -22.77   <2e-16 ***
## result$FCH4  4.567e-01  3.026e-03  150.94   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006543 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.6999, Adjusted R-squared:  0.6999 
## F-statistic: 2.278e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.049401 -0.003553  0.000473  0.003228  0.215192 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0024528  0.0001712  -14.33   <2e-16 ***
## test_set$FCH4  0.1162034  0.0073033   15.91   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008413 on 2441 degrees of freedom
## Multiple R-squared:  0.09397,    Adjusted R-squared:  0.0936 
## F-statistic: 253.2 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.193797 -0.002251  0.000303  0.002044  0.199596 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.503e-03  7.452e-05  -20.18   <2e-16 ***
## result$FCH4  4.641e-01  3.382e-03  137.24   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007313 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.6585, Adjusted R-squared:  0.6584 
## F-statistic: 1.884e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.060426 -0.003351  0.000501  0.003240  0.050792 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0025748  0.0001328  -19.39   <2e-16 ***
## test_set$FCH4  0.1206494  0.0060113   20.07   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006524 on 2441 degrees of freedom
## Multiple R-squared:  0.1416, Adjusted R-squared:  0.1413 
## F-statistic: 402.8 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.112170 -0.002196  0.000342  0.002042  0.076036 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.609e-03  6.172e-05  -26.07   <2e-16 ***
## result$FCH4  4.464e-01  2.801e-03  159.36   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006057 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.7222, Adjusted R-squared:  0.7222 
## F-statistic: 2.54e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.042367 -0.003556  0.000441  0.003464  0.058800 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0027346  0.0001368  -19.99   <2e-16 ***
## test_set$FCH4  0.0801529  0.0058239   13.76   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00671 on 2441 degrees of freedom
## Multiple R-squared:  0.07201,    Adjusted R-squared:  0.07163 
## F-statistic: 189.4 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.202201 -0.002323  0.000209  0.002042  0.193362 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.480e-03  7.423e-05  -19.94   <2e-16 ***
## result$FCH4  4.638e-01  3.369e-03  137.69   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007285 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.6599, Adjusted R-squared:  0.6599 
## F-statistic: 1.896e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.068185 -0.003803  0.000247  0.003104  0.298659 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0023083  0.0002276  -10.14   <2e-16 ***
## test_set$FCH4  0.1388577  0.0100618   13.80   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01118 on 2441 degrees of freedom
## Multiple R-squared:  0.07238,    Adjusted R-squared:  0.072 
## F-statistic: 190.5 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.163160 -0.002265  0.000230  0.001939  0.288026 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.442e-03  7.978e-05  -18.08   <2e-16 ***
## result$FCH4  4.877e-01  3.620e-03  134.71   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007829 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.6501, Adjusted R-squared:   0.65 
## F-statistic: 1.815e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.048476 -0.003300  0.000522  0.003111  0.119566 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0025336  0.0001343  -18.86   <2e-16 ***
## test_set$FCH4  0.1209069  0.0064598   18.72   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006589 on 2441 degrees of freedom
## Multiple R-squared:  0.1255, Adjusted R-squared:  0.1251 
## F-statistic: 350.3 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.173352 -0.002135  0.000328  0.001966  0.094857 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.514e-03  6.244e-05  -24.24   <2e-16 ***
## result$FCH4  4.699e-01  2.834e-03  165.84   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006128 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.7379, Adjusted R-squared:  0.7379 
## F-statistic: 2.75e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.043127 -0.003158  0.000817  0.003417  0.036764 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0029084  0.0001279  -22.75   <2e-16 ***
## test_set$FCH4  0.0929055  0.0048287   19.24   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006307 on 2441 degrees of freedom
## Multiple R-squared:  0.1317, Adjusted R-squared:  0.1313 
## F-statistic: 370.2 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.193853 -0.002244  0.000617  0.002468  0.076347 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.979e-03  6.651e-05  -29.75   <2e-16 ***
## result$FCH4  3.900e-01  3.018e-03  129.22   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006527 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.6309, Adjusted R-squared:  0.6308 
## F-statistic: 1.67e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.063779 -0.003421  0.000621  0.003434  0.276040 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0025469  0.0001862 -13.676   <2e-16 ***
## test_set$FCH4  0.0763056  0.0086258   8.846   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.009131 on 2441 degrees of freedom
## Multiple R-squared:  0.03106,    Adjusted R-squared:  0.03067 
## F-statistic: 78.25 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.140648 -0.002291  0.000272  0.001972  0.263238 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.391e-03  7.863e-05   -17.7   <2e-16 ***
## result$FCH4  4.922e-01  3.568e-03   137.9   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007716 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.6608, Adjusted R-squared:  0.6607 
## F-statistic: 1.903e+04 on 1 and 9769 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.063344 -0.003121  0.000654  0.003023  0.117039 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -0.0024969  0.0001436  -17.39   <2e-16 ***
## test_set$FCH4  0.1065710  0.0061119   17.44   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007044 on 2441 degrees of freedom
## Multiple R-squared:  0.1108, Adjusted R-squared:  0.1104 
## F-statistic:   304 on 1 and 2441 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.182305 -0.002217  0.000348  0.001957  0.087063 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -1.511e-03  6.643e-05  -22.75   <2e-16 ***
## result$FCH4  4.379e-01  3.015e-03  145.25   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.006519 on 9769 degrees of freedom
##   (19070 observations deleted due to missingness)
## Multiple R-squared:  0.6835, Adjusted R-squared:  0.6835 
## F-statistic: 2.11e+04 on 1 and 9769 DF,  p-value: < 2.2e-16